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We have performed 31 P-NMR measurements on Ce(Rui_ x Fe cc )PO in order to investigate ferro- 
magnetic (FM) quantum criticality, since heavy-fermion (HF) ferromagnet CeRuPO with a two- 
dimensional structure turns to a HF paramagnet by an isovalent Fe- substitution for Ru. We found 
that Ce(Ruo.i5Feo.85)PO shows critical fluctuations down to ~ 0.3 K, as well as the continuous 
suppression of Curie temperature and the ordered moments by the Fe-substitution. These exper- 
imental results suggest the presence of a FM quantum critical point (QCP) at x ~ 0.86, which is 
a rare example among itinerant ferromagnets. In addition, we point out that the critical behaviors 
in Ce(Ruo.i5Feo.85)PO share the similarity with those in YbRli2Si2, where the local criticality of 
/ electrons has been discussed. We reveal that Ce(Rui_ x Fe x )PO is a new system to study FM 
quantum criticality in HF compounds. 

PACS numbers: 76.60.-k, 75.30.Kz 74.70.Xa 



It has been believed that a ferromagnetic (FM)- para- 
magnetic (PM) transition at zero magnetic field is a text- 
book example of a 2nd-order phase transition at finite 
temperatures However, there have been a lot of re- 
ports that the FM-PM transition changes from 2nd-order 
to lst-order at a tricritical point (TCP) in metallic fer- 
romagnets, when the FM transition is suppressed to K 
by non-thermal control parameters [2 -6]. In addition, the 
quantum FM transition was predicted to be generically of 
lst-order from the theoretical viewpoints 0, HI- Nature of 
the quantum phase transition from the FM to PM state 
in metallic ferromagnets has been revisited recently. 

A quantum critical point (QCP), defined as a zero- 
temperature 2nd-order phase transition driven by non- 
thermal physical parameters, has attracted much atten- 
tion, since fascinati ng p henomena, e.g., unconventional 
superconductivity 0,[lQj and nematic order [TlJ have been 
discovered near a QCP. For studying the nature of a QCP, 
heavy-fermion (HF) systems, in which / electrons behave 
as itinerant electrons with heavy effective-mass, are one 
of the best systems for this purpose, since a magnetic 
transition can be tuned by various physical parameters 
such as chemical substitution H, [l3j , pressure 3 [Hj or 
mag netic field[16]. Up to now, there have been consider- 
able efforts for understanding the nature of a QCP on var- 
ious antiferromagnetic (AFM) HF compounds, but very 
few on FM HF compounds, particularly on Ce-based FM 
HF compounds. 

The iron oxypnictide CeFe(Ru)PO is a related material 
of the iron-based superconductor LaFePO[l7j- They pos- 
sess the same two-dimensional layered structure, stacking 
the Ce(La)0 and Fe(Ru)P layers alternatively. The CeO 
layer contributes to the large magnetic response, and the 
Fe(Ru)P layer is conductive in CeFe(Ru)PO. CeRuPO 
is a FM HF compound with Curie temperature Tq = 
15 K and coherent temperature 7k — 10 K [18]. On 



the other hand, counterpart CeFePO is a HF compound 
with a PM ground state down to 80 mK and possesses a 
large Sommerfeld coefficient 7 = 700 mJ/mol-K 2 at low 
temperatures [19]. Quite recently, we reported that Ce- 
FePO shows a metamagnetic (MM) behavior related to 
the Kondo effect when magnetic fields are applied per- 
pendicular to the c axis, and exhibits a non Fermi liquid 
behavior at the MM field jioj . Therefore, we expect that 
the ground state can be tuned continuously from FM to 
PM by substituting isovalent Fe for Ru. Actually, con- 
tinuous change from FM to PM state revealed by mag- 
netization and specific heat measurements down to 2 K 
is reported in CeFe(Asi_ :E P :E )0 [ll|. In particular, we 
anticipate that Ce(Rui_ :c Fe ;E )PO is an ideal system for 
investigating physical properties of FM quantum critical- 
ity from a microscopic point of view by 31 P-NMR mea- 
surements, since the substitution is done at the Fe(Ru)P 
layer and thus introduces small randomness effect at the 
magnetic CeO layer as well as at the P site detected with 
NMR. Here we report evolution of magnetic properties in 
Ce(Rui_ ;E Fe ;r )PO revealed by 31 P-NMR measurements 
and show the phase diagram in this system. We stress 
that an NMR measurement is one of the best experi- 
mental technique for distinguishing between 1st and 2nd 
order transition as mentioned in our previous paper [22] . 

The polycrystalline Ce(Ru i _ :c Fe :c )PO (x = 0.25, 0.5, 
0.75, 0.85, 0.87, 0.9, 0.95, and 1) was synthesized by 
solid-state reaction [H, Hif . As reported in the previous 
paper [2o|, these possess a two-dimensional anisotropy, so 
the polycrystalline samples were uniaxially aligned by 
taking advantage of the anisotropy of the magnetic sus- 
ceptibility. The c axis of the samples are mostly aligned, 
but the a and b axes are randomly oriented in the ab 
plane. 31 P-NMR measurements were performed on the 
aligned samples. 

NMR spectra were obtained with the field-swept tech- 
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FIG. 1: (Color online) (a) T dependence of Mi/H evaluated 
from the Knight shift determined at the peaks of H J_ c 
and H \\ c spectra obtained at 0.5 T for Ce(Rui_ccFe x )PO 
(x = 0.25, 0.5, 0.75, 0.85, 0.87, 0.9, 0.95, and 1). The Knight 
shift K± and K c are also represented in the right axis. Mi/H 
is estimated from the equation Mi/H = (Ki — i^ rb)/^.hf,iso 
where isotropic hyperfine coupling constant 31 Ahf,i so = 0.2 
T/fiB and T-independent orbital part of Knight shift K OT h 
= 0.09 % derived at x = 1 were used. Below x = 0.85, Ki 
increases below Tc due to the appearance of the internal field 
at the P site. The arrows indicate Tc. (Inset) T dependence 
of Mi/H above x = 0.87. M±/H shows a peak at around 
T max ~ 5 K (indicating the arrows) instead of the FM order, 
suggestive of the PM HF ground state at low fields, (b) T 
dependence of Mi/H (upper panel) and in-plane dynamical 
susceptibility S± probed by 1/Ti (lower panel) at x — 0.25. 
Tc (broken line) is unambiguously determined by the peak of 
S±. 

nique, and the Knight shift, which is the measure of the 
local susceptibility at the nuclear site, was determined 
from the peak field of each resonance spectrum. Knight 
shift Ki(T, H) (i =_L and c) is defined as, 
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where H res is a magnetic field at a resonance peak, and 
Hq and ujo are the resonance field and frequency of a bare 
31 P nucleus and have the relation of ujq = 7 n i^o with the 
31 P-nuclear gyromagnetic ratio j n . The Knight shift was 
measured in H parallel (K c ) and perpendicular (K±) to 
the c axis for all samples at ji^H ~ 0.5 T. Reflecting the 
two-dimensional XF-type anisotropy, K c proportional to 
the out-of-plane component of the susceptibility is small 
and almost T independent, but K± proportional to the 
in-plane susceptibility shows strong T dependence origi- 
nating from the Curie- Weiss (CW) behavior of x(T) at 
high temperatures as shown in Fig.[]Ja). At x = 0.25, K c 
suddenly increases below Tq = 10 K due to the appear- 
ance of the internal field at the P site, but K± gradually 
increases from much higher temperature than Tc due to 
the CW behavior. Tq is unambiguously determined by 
the in-plane low-energy spin fluctuations S± as shown in 
Fig. DJb) probed by the nuclear spin-lattice relaxation 
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FIG. 2: (Color online) T dependence of low-energy in-plane 
spin fluctuations S± at ~ 0.5 T. (See in the text.) Below 
x = 0.85, T dependence of S± shows a peak at Tc (indicating 
the arrows). Above x = 0.9, S± as well as K± along both 
directions becomes almost a constant at low temperatures, 
indicative of the formation of a HF ground state, although 
S± continue to increase down to 0.08 K at x = 0.87. (Inset) 
x dependence of S± at 0.2 K. S± diverges toward x ~ 0.86. 
A broken curve is guide to eyes. 

rate (1/Ti ) in H || c. As discussed in the previous pa- 
per [20[, 1/Ti probes spin fluctuations perpendicular to 
the applied magnetic field, and thus 1/Ti in H \\ c and 
H _L c are described as, 



f) -My 



= S C + S ± . 



(2) 
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From two equations, we can know low-energy g-summed 
spin fluctuation of in-plane (S±) and out-of-plane (S c ) 
component, separately. 

By substitution of Fe for Ru, Tq and the internal field 
at the P site, proportional to ordered magnetic moment 
(/x rd)j continuously approach to zero toward x ~ 0.86, 
suggesting that a FM QCP is located at x ~ 0.86. The 
x dependence of S± also indicates the existence of a FM 
QCP at x - 0.86 as shown in Fig. [2 

Above x = 0.9, spin fluctuation S± as well as K± be- 
comes almost constant at low temperatures, indicative of 
the formation of the PM HF state, although S± continue 
to increase down to 0.08 K at x = 0.87. At x = 1, as 
mentioned above, we reported MM behavior at the MM 
field Hm — 4.3 T, and suggested that Hm corresponds 
to an energy breaking the HF state, since the linear re- 
lationship holds in various Ce-based HF metamagnets, 
between Hm and the temperature where susceptibility 
shows a broad maximum T max or the inverse of the effec- 
tive mass. We investigated the x revolution of the MM 
behavior below x = 0.85. Figure [31(a) shows the H de- 
pendence of magnetization M± at 0.1 K below x = 0.87 
and at 0.5 K for x = 0.85 just above Tc, which are PM 
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FIG. 3: (Color online) (a) H dependence of magnetiza- 
tion Mj_ (if) in H ± c at 0.1 K below x = 0.87 and at 
0.5 K for x = 0.85 just above Tc using the relation of 
M±(H) = K±(H)H/ Am, iso- Dashed curves are guide to eyes. 
Above x = 0.87, the M±(H) suddenly increases with increas- 
ing H and deviates from linear relation in the field range of 
1 - 5 T, which is a definition of a metamagnetic behavior, 
although M±(H) increase linearly at x = 0.85. (b) H de- 
pendence of (1/TiT)h±c at 0.1 K below x = 0.87 and at 
0.5 K for x = 0.85. At x = 1, the low-field PM state with 
HF character is separated from the high-field polarized PM 
state by a broad peak of (1/TiT)h±c at Hm — 4.3 T. Sim- 
ilar broad peak (1/TiT)h±c was observed at x = 0.95, but 
is changed to a shoulder structure at x = 0.9. On the other 
hand, (1/TiT)h±c is monotonically suppressed by magnetic 
field and no anomaly is observed at x = 0.85 and 0.87. We 
define Hm from the H dependence of M± and (1/TiT)h±c as 
shown by the dotted lines. Hm continuously decreases with 
decreasing x. Dashed curves are guide to eyes. 

state at zero field. M± is defined as M± = K±H/A^^ ao 
with the isotropic hyperfine coupling constant 31 Ahf,i so 
= 0.2 T//iB derived at x = l[20j. We observed the sim- 
ilar MM behavior at x = 0.87,0.9 and 0.95, and found 
that Hm decreases with decreasing x, consistent with the 
suppression of T max with decreasing x. [See in the inset 
of Fig. HJa).] The MM behavior was also confirmed from 
the spin dynamics. Figure [3fb) shows the H dependence 
of (1/TiT)h± c - As reported previously [20|, at x = 1, 
the low-field PM state with the HF character is sepa- 
rated from the high-field polarized PM state with small 
(l/T 1 T) H ±c by a broad peak of (1/T 1 T) H ± C at H M ~ 
4.3 T. Similar broad peak (1/TiT)h±c was observed at 
x = 0.95, but is changed to the shoulder structure at 
x = 0.9. On the other hand, (1/TiT)h± c is monoton- 
ically suppressed by magnetic field and no anomaly is 
observed at x = 0.85 and 0.87. We determined Hm from 
the H dependence of M± and (1/TiT)h±c as shown by 
the dotted lines in Fig. 02 At x = 0.85, we could not 
observe any MM behavior both in M± and (l/T 1 T) H ±c 



FIG. 4: (Color online) x — T — H phase diagram in 
Ce(Rui_ x Fe ;E )PO determined by NMR results. T c in 
CeRuPO is continuously suppressed by Fe substitution, and 
approaches to zero at x ~ 0.86, suggestive of the presence of 
a FM QCP at x ~ 0.86. Furthermore, Hm observed in Ce- 
FePO decreases with decreasing x and seems to merge with 
the FM QCP at x ~ 0.86. (Inset) x dependence of ordered 
magnetic moment (/i or d) estimated from K c . (/i or d) continu- 
ously approaches to zero toward x ~ 0.86, which is consistent 
with x dependence of Tc- Dashed curve is guide to eyes. The 
2nd-order FM transition is suggested from the development 
of the critical fluctuations down to low temperatures as well 
as the continuous suppression of Tc and the ordered moments 

(Mord)- 

On the basis of above results, we developed the x — T — 
H phase diagram on Ce(Rui_ :E Fe :c )PO in Fig. [4j Tq in 
CeRuPO is continuously suppressed by the Fe substitu- 
tion, and approaches to zero at x ~ 0.86. This is consis- 
tent with the x dependence of (/x rd) evaluated from the 
K C1 which is shown in the inset of Fig. 2J In addition, 
remarkable divergence of S± against T at x = 0.85 and 
against x at 0.2 K as shown in Fig. [2] strongly suggest the 
presence of a FM QCP at x - 0.86 in Ce(Ru i _ ;r Fe :c )PO. 
Furthermore, Hm observed in CeFePO decreases with 
decreasing x and seems to merge with the FM QCP at 
x ~ 0.86, suggestive of the close relationship between the 
FM QCP and the criticality of the MM behavior. 

There have been a lot of theoretical and experimen- 
tal efforts for understanding a FM quantum critical- 
ity, since a lst-order FM transition has been observed 
in various itinerant ferromagnets, e.g. ZrZn2, MnSi, 
(Sri_ x Ca :c )Ru03 and so on[25[. In addition, Belitz et al. 
proposed that the lst-order MM transitions emerge from 
the TCP in the small field range and terminate at quan- 
tum critical end points [§[. This characteristicphase dia- 
gram was experimentally confirmed in UGe2 [5] , URhGe 
Q and partly in UCoAl(26|. The phase diagram ob- 
tained in Ce(Rui_ a: Fe :E )PO is quite different from that in 
other itinerant ferromagnets reported so far [e.g. UGe2 
and Co(Si_a;Se x )2@|]. Particularly, the 2nd-order char- 
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FIG. 5: (Color online) T dependence of M ± /H (K±) (a) 
and S± (b) measured at the smallest field we could ( fioH ~ 
0.07 T) at x = 0.85. Static susceptibility related with q = 
becomes saturate below 3 K, although S± continues to de- 
velop toward Xc , suggesting that spin fluctuations would pos- 
sess AFM (q y£ 0) components as well as the FM (q = 0) com- 
ponents near the FM QCP. These behaviors were in sharp con- 
trast with those observed in x < 0.75 as shown in Fig. [T](b), 
but were actually observed in YbRli2Si2 near an if-induced 
AFM QCP, in which the possibility of the local criticality has 
been discussed. 



acter of the FM QCP is a rare example among itinerant 
ferromagnets [27) . However, it might be possible to rec- 
oncile the different phase diagram by the suppression of 
a TCP temperature to nearly zero or a negative value 
in Ce(Rui_ :E Fe :E )PO. Absence of a lst-order MM tran- 
sition is also consistent with the 2nd-order character of 
the FM-PM transition. 

In addition to the novel 2nd-order FM QCP, it de- 
serves to be noted that the "wing" of the MM crossover 
also converges at x ~ 0.86. Since the Hm is related to 
the characteristic energy of the HF state, this reminds 
us of the local QCP observed in YbRh 2 Si 2 [i!, Ac- 
tually, similarity was observed between critical behaviors 
in YbRh2Si2 and those in Ce(Ruo.i5Feo.85)PO. Figure [5] 
shows the T dependence of M±/H (K±) and S± mea- 
sured at the smallest field we could ( ji^H ~ 0.07 T) 
at x = 0.85. Static susceptibility related with q = 
becomes saturate below 3 K, although S± continues to 
develop with a S± oc T~ l relation toward Xc, suggesting 
that spin fluctuations would possess AFM (q ^ 0) com- 
ponents as well as the FM (q = 0) components near the 
FM QCP. These behaviors were in sharp contrast with 
those observed in x < 0.75 [see in Fig. OJb)], but were 
actually observed in YbRh 2 Si 2 near AFM QCP0. It is 
noted that both compounds possess the two-dimensional 
XY anisotropy in the PM state, which might be related 
to the origin of the local QCP. For the local QCP, a sud- 



den change of the Fermi surface is expected at the QCP, 
where magnetic fluctuations at various q components are 
considered to develop due to the Fermi-surface instabil- 
ity. Thus, we speculate that sudden Fermi-surface change 
might occur at x ~ 0.86, which cannot be understood by 
the ordinary spin density wave scenario [li l3lL l32| . 

In summary, both of the static and dynamic suscep- 
tibilities in the FM and PM states of Ce(Rui_ ;r Fe ;r )PO 
clearly suggest that the FM QCP is located at x ~ 0.86, 
and that the FM QCP possesses a 2nd-order character, 
which is a rare example for itinerant ferromagnets. In 
addition, continuous suppression of the MM crossover 
toward x ~ 0.86, and the characteristic behavior of S± 
suggests that the close relationship between the FM QCP 
and the criticality of the Kondo effect, which can be un- 
derstood by the local QCP scenario. 
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